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Abstract Compared with other terrestrial environments, 
the stream environment generally presents a linear spatial 
structure and relatively simple environment. In a stream 
landscape, the dispersal direction of stream-type organisms 
usually presents a linear structure along the stream, 
which results in the limited dispersal and the genetic 
differentiation of stream-type organisms across different 
stream sections. The Shangcheng Stout Salamander 
(Pachyhynobius shangchengensis) is a narrowly distributed 
stream salamander in Dabie Mountains of East China. 
In the present study, we tested for the impact of stream 
landscape (i.e. waterfalls and underground river) on 
genetic structure and dispersal pattern in P. shangchengensis 
based on 12 nuclear microsatellite loci from 195 individuals 
in 3 populations (A, B and C) from three closely connected 
sections within one stream. Genetic diversity results 
suggested that Population B contains relatively high 
genetic diversity for P. shangchengensis when compared 
to the other populations (A and C). Detectable genetic 
differentiation was found (Իլ = 0.008, P = 0.007) among 
three populations, which was also supported by the 
Structure, FCA analysis and relatedness estimates of each 
pair of individuals among populations. The assignment 
test suggested that P. shangchengensis has philopatric males 
and female-biased dispersal (mean female Alc = -0031, SE 
= 0.225; mean male Alc = 0026, SE = 0.198). Female-biased 
dispersal was also supported by analyses for each sex 
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(ie. Spatial autocorrelation, Genetic distance, Relatedness 
analysis). Our study indicated that small and isolated 
populations (A and C) had relatively low genetic diversity 
due to the limited population size. For stream salamanders, 
landscape features (16. waterfalls and underground river) 
can influence the ability of an individual to disperse through 
the landscape, and consequently influence the formation of 
strong genetic differentiation of P. shangchengensis. 


Keywords Dabie Mountains, genetic structure, 
landscape features, Pachyhynobius shangchengensis, sex- 


biased dispersal 


1. Introduction 


Streams are fractal-like, occur with a hierarchical linear 
structure where smaller stream channels join to form 
larger ones in a dendritic pattern that resembles branches 
on a tree (Horton, 1945). It is because of the hierarchical 
linear structure of streams that environments and 
resources differ across the dendritic streams. The spatial 
heterogeneity of this environment and resources have led 
to discrete transitions in morphological, physiological, 
and behavioral traits (Wilbur, 1980; Moyle and Cech, 
1988; Duellman and Trueb, 1994; Lowe et al., 2004; 
Mullen et al., 2010). It may have facilitated the formation 
of different dispersal patterns and genetic differentiation 
within different stream sections (Finn et al., 2006; 
Muneepeerakul et al., 2008; Grant et al., 2009; Mullen 
et al., 2010). Therefore, in a stream landscape, studies 
focusing on how factors (e.g. resource richness, stream 
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landscape characteristics) have affected landscape-scale 
population structure and scaling of dispersal are critical 
for testing general models (Lowe et al., 2006; Campbell 
Grant et al., 2007; Mullen et al., 2010). 

Intraspecific interactions, competitive ability, habitat 
characteristics, and resource distribution and availability 
can influence dispersal patterns (Clobert et al., 2001; 
Handley and Perrin, 2007). The dispersal patterns of 
species function as the key influence on genetic structure 
and help maintains genetic diversity (Bowler and Benton, 
2005; Clobert et al., 2001; Orsini et al., 2013; Smith et al., 
2016; Olah et al., 2017). When the evolutionary pressures 
between sexes differ, asymmetric dispersal of the sexes 
may occur, called sex-biased dispersal (Bonte et al., 2011). 
Several theories have been generated to explain the 
differences in dispersal patterns between different sexes, 
such as local resource competition theory (Greenwood, 
1980), local mating competition theory (Dobson, 1982) 
and avoidance of inbreeding theory (Pusey, 1987). 
In vertebrates, mammals usually show male-biased 
dispersal and birds typically show female-biased dispersal 
(Greenwood, 1980; Pusey, 1987; Handley and Perrin, 
2007). As to the sex-biased dispersal reviewed in fish 
(Bekkevold et al., 2004; Stiver et al., 2007), amphibians 
(Palo et al., 2004; Liebgold et al., 2011; Wang et al., 2012) 
and reptiles (Dubey et al., 2008), the dispersal patterns and 
mating strategies within these taxa are diverse. 

Unlike other landscape types, in stream landscapes, 
stream-type organisms are usually confined to rivers and 
streams, which leads to present a linear dispersal pattern 
along the waterways. As described in the above paragraph, 
previous studies have shown that some stream-dwelling 
organisms own different evolutionary, demographic and 
ecological processes across the dendritic streams (Finn et 
al., 2006; Muneepeerakul et al., 2008; Grant et al., 2009; 
Mullen et al., 2010). Therefore, to provide important 
evidence for the correlation between the stream landscape 
and the formation of different evolutionary, demographic 
and ecological processes in hierarchical linear structure, 
more work is needed that focuses on the dispersal 
patterns of stream-type organisms. 

Population connectivity is critical to many ecological 
processes and conservation goals, and can be a key 
factor affecting the regional viability of populations 
(Curtis and Bevers, 2002; Hayden, 2006; Saastamoinen 
et al., 2018). For amphibians, there are several types of 
landscape features, such as the presence of wetlands, 
changes in slope and cover type, that should influence 
the dispersal and population genetic structure (Spear 
et al., 2005). Amphibians generally require wet areas to 
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survive and reproduce (Duellman and Trueb, 1994), and 
therefore the presence of wetlands should be important 
for maintaining gene flow. Differences in slope among 
sites can account for the changes in dispersal or gene flow, 
as high levels of amphibian population subdivision exist 
among populations separated by mountains compared to 
low subdivision within basins in these populations (Funk 
et al., 2005), ostensibly due to increased slope and reduced 
mobility. In fact, the effects of slopes and topography 
have been identified as the main predictors of amphibian 
dispersal or gene flow (Funk et al., 2005; Olah et al., 2017). 
Cover type may have a strong influence on landscape 
connectivity for amphibians. Intermittent landscape 
features directly affect population dispersal or gene 
flow (Spear et al., 2005; Talbot et al., 2017). Landscape 
genetics uses genetic data as the dependent variable to 
correlate genetic relationships with several independent 
variables representing natural, physical obstacles to 
address some of the most challenging ecological and 
evolutionary questions (e.g., source-sink dynamics, 
barriers and corridors, influence of landscape change) 
(Storfer et al., 2007; Saastamoinen et al., 2018). In a linear 
stream landscape, there are many landscape factors, such 
as waterfall, underground river, changes in slope and 
cover type, that can cause the stream discontinuity. For 
stream-type organisms, linear hierarchical networks 
with landscape barriers may prevent local adaptive 
populations from being submerged and promote genetic 
differentiation (Lowe et al., 2006; Campbell Grant et al., 
2007; Campbell Grant et al., 2010; Mullen et al., 2010). 
The Shangcheng Stout Salamander (Pachyhynobibus 
shangchengensis) is a narrowly distributed salamander 
endemic to streams above 500 meters in the Dabie 
Mountains, Eastern China (Fei et al., 2010). It has deep 
genetic divergence among regional populations on a 
small spatial scale (Zhao et al., 2013; Pan et al., 2014, 
2019). For narrowly distributed species with deep genetic 
divergence, such as P. shangchengensis, how it dispersal 
across those nature physical obstacles, such as waterfalls 
and underground rivers, in a linear river system, will 
help the understanding of how physical obstacles allow 
or restrict dispersal and create opportunities for genetic 
differentiation, then led to the deep genetic divergence. 
Therefore, in this study, we use an individual-based 
genetic approach with 12 microsatellite loci to assess how 
stream landscape features (e.g. waterfall, underground 
river) potentially influence the dispersal pattern and 
genetic structure of P. shangchengensis in Yaoluoping 


National Nature Reserve of Dabie Mountains. 
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2. Materials and Methods 
2.1. Ethics statements The collection of samples was 
performed within a long-term investigation project on 
amphibians of Dabie Mountains. This investigation 
project and the sample collection were approved by Anhui 
Yaoluoping National Nature Reserve, Anhui Province, 
China. Additionally, the protocol was approved by the 
ethics committee of Anhui Normal University. 


2.2. Study area and sampling Our study was conducted 
in the Yaoluoping National Nature Reserve (30°57'57" 
N, 116°4'5" E) in Southern Dabie Mountains, Yuexi 
county, Anhui Province in Eastern China (Figure 
1). Dabie Mountains are located in the ecotone of 
subtropical evergreen broad-leaved forest and warm- 
temperate deciduous broad-leaved forest zone (Zheng 
et al., 2012). Pachyhynobius shangchengensis is widely 
distributed throughout the streams and has a relatively 
large population in Yaoluoping National Nature Reserve, 
which will facilitate the population recovery after our 
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sampling. 

From August 26 to August 28, 2013, our study was 
carried out in three river sections (about 600 m; A, B 
and C) of the Yangtianwo valley area in the Yaoluoping 
National Nature Reserve (Figure 2). The streams 
have slight alkalinity (pH of 7.3), moderate midday 
temperatures (15.7 °C), river flow velocity (0.028-0.669 m/ 
s). The B and C sections are part of the main river, which 
are isolated from each other by a waterfall (15 m height); 
section A is a tributary of the main river. The landscape 
between sections A and C is an underground river, which 
is mostly covered by vegetation section (10 m) without a 
pool of water. The landscape details for sampling along 
the river are shown in Figure 2. 

All adult individuals were captured by dip netting with 
relative positions recorded with a 500-meter Engineering 
ruler. Muscle or liver tissue was sampled and preserved in 
100% ethanol for DNA extraction later. Specimens were 
fixed in 10% formalin and subsequently transferred to 


75% ethanol for storage. All the specimens were deposited 


116°10'0"E 


Figure 1 Sampling area for the populations of Pachyhynobibus shangchengensis. The red line represents the county line. Values with different colors rep- 
resent altitudes. The green star corresponds to sample areas. Huoshan County and Yuexi County are located in the Anhui Province. Yingshan County is 


located in the Hubei Province. 
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Figure 2 The schematic diagram of the sampled river in Yang- 
tianwo of Yaoluoping area, Yuexi Country, China. The width of 
the sampled river is no more than 2 meters. The B and C sections 
are part of the main river, which are isolated from each other by 
a waterfall (15 m height); section A is the tributary of main river. 
The landscape between sections A and C is the underground river 
(10 m). 
at the Zoological Specimen of Museum Anhui University 
(AHU), Hefei, Anhui, China. Sex was determined through 
sexually dimorphic features between sexes, such as the 
head length and width (Clemen and Greven, 2009) and the 
presence of reproductive system (Bulakhova and Berman, 


2013; Vadim et al., 2015). 


2.3. Lab work Total DNA was extracted from samples 
using a standard proteinase K/ phenol-chloroform 
protocol (Sambrook et al., 1989). An EasyPure PCR 


Purification Kit (TransGene) was used to purify each 
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DNA extraction. Twelve independent polymorphic 
microsatellite loci for P. shangchengensis were selected 
from previous studies (Table 1) (Wang et al., 2010). PCR 
was performed in 20 pL reactions containing 20-50 ng 
of template DNA, 10 uL 2 x EasyTaq PCR SuperMix 
polymerase (TransGen Biotech, containing 1.25U Ex 
Taq, 0.4 mmol/L dNTP, 4 mmol/L Mg”) and 0.2 pmol/L 
of each primer. Thermal cycling was performed with the 
following profile: 5 min denaturing step at 95 °C; followed 
by 35 cycles of 30s at 95 °C, 1 min at optimal annealing 
temperature (Wang et al., 2010) and 1 min at 72 °C; and 
a final extension at 72 °C for 10 min. All PCR products 
were visualized on an ABI 3730 sequencer (PE Applied 
Biosystems) with GS500 marker. Allele peaks were scored 
using GeneMarker 1.85 (SoftGenetics LLC) (Holland and 
Parson, 2011). 


2.4. Basic genetic information Number of alleles (Na), 
effective number of alleles (Ne), observed heterozygosity 
(Ho), expected heterozygosity (H,) and fixation index (F) 
were estimated in GenAIEx 6.5 (Peakall and Smouse, 
2006). Cervus 3.0 was used to assess the polymorphism 
information content (PIC) and the frequency of null alleles 
(Kalinowski et al., 2007). 


2.5. Population genetic structure To determine 
whether there is genetic substructure among these 
three populations (A, B and C), we performed a genetic 
clustering analysis using STRUCTURE v.2.3 (Pritchard 
et al., 2000). The lengths of Markov Chain Monte Carlo 
iterations and burn-in were set as 1,000,000 and 100,000, 
respectively. We tested the range of possible clusters (K) 
tested from 1 to 10, and performed 10 independent runs 
for each K. The other parameters were set as following: 
admixture model and correlated allele frequencies 
between populations. 

Analysis of molecular variance (AMOVA) in GenAlEx 
was used to partition genetic variation within and among 
these sampling sites and to estimate overall and pairwise 
population genetic differentiation (Fr) (Wright, 1965; 
Excoffier et al., 1992; Peakall and Beattie, 1995). In order 
to assess if there was any bias in Fy; estimation given the 
uneven sample size in Population B, we performed an 
AMOVA on randomly subsampled sets of 20 samples from 
each population using the ‘shuffle’ function of GenAlEx. 
Tests for genetic differentiation were performed by 1000 
random permutations. For probing the genetic structure, 
we conducted the factorial correspondence analysis (FCA) 
by the GENETIX 4.05 (Belkhir et al., 1996-2004). The 
genetic differentiation among geographical samples was 
calculated by the standardized statistics Ի (Meirmans, 
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Table 1 Average genetic variation of twelve polymorphic microsatellite loci for Pachyhynobibus shangchengensis. 


Locus Size range Na Ne Ho 
PSH29 160-193 15 7.49 0.84 
PSH74 167-271 21 13.19 0.86 
PSH130 229-273 13 7.41 0.68 
PSH236 147-207 16 5.56 0.77 
PSH133 199-291 18 421 0.49 
PSH255 97-161 17 10.56 0.85 
PSH422 126-204 19 11.29 0.81 
PSH463 168-228 15 10.86 0.73 
PSH493 156-212 13 8.40 0.81 
PSH146 184-270 27 13.02 0.90 
PSH118 161-177 10 4.41 0.48 
PSH131 116-138 6 1.58 0.29 
Mean 15.83 8.17 0.71 


209 

H; F Fs PIC HWE 
0.87 0.03 0.07 0.85 

0.92 0.07 0.10 0.92 

0.87 0.21 0.14 0.85 “A 
0.82 0.06 0.09 0.80 

0.76 0.36 0.30 0.73 iG 
0.91 0.06 0.001 0.89 

0.91 0.12 0.19 0.91 

0.91 0.20 0.13 0.90 “A 
0.88 0.08 0.02 0.86 

0.92 0.02 0.03 0.92 

0.77 0.38 0.46 0.73 ਨ 
0.37 0.20 0.23 0.33 “A 
0.83 0.15 0.15 0.81 


Notes: Presented are locus code, fragment size ranges, number of alleles (Na), effective number of alleles (Ne), observed heterozygosity (Ho), expected 
heterozygosity (H,) and fixation index (F). The populations significant departure (P Հ 0.01) from HWE. 


2006) and Jost's Dest (Meirmans and Hedrick, 2011) 
included in GenAlEx. 


2.6. Sex-biased dispersl analyses Assignment tests 
implemented in GenAlEx (Favre et al., 1997; Mossman 
and Waser, 1999) were used to investigate the pattern of 
sex-biased dispersal in P. shangchengensis. A log likelihood 
assignment test value is calculated for each individual 
and then the Assignment Index correction (AIC) score 
is calculated by subtracting the mean log likelihood of 
the population from an individual’s log likelihood score. 
For each population Alc values will average zero, while 
individuals with negative Alc scores are considered to 
have a higher probability of being immigrants (Mossman 
and Waser, 1999). The significance of sex-biased dispersal 
was examined using Mann Whitney U tests according to 
the difference between sexes in the frequency distribution 
of Alc scores. A negative score indicates male-biased 
dispersal and a positive value is indicative of female- 
biased dispersal (Mossman and Waser, 1999). The females 
should be more related to proximate females than distant 
females when females are philopatric, and the opposite 
pattern should occur if males are philopatric. 

To further evaluate the inbreeding pattern within 
or among populations, we calculated the relatedness 
estimate (r) (Wang, 2002) and the inbreeding coefficient 
(Fis) for each pair of individuals (Ritland, 1996) by 
COANCESTRY 1.0.1.1 (Wang, 2011). In addition, we also 
estimated the pairwise relatedness values (Queller and 
Goodnight, 1989) and pairwise relationship (Kalinowski 
et al., 2006) based on the allele frequencies for determining 
the distribution pattern of kinship among individuals 


within Population B by KINGROUP v.2 (Konovalov et 
al., 2004). In addition, for further probing spatial genetic 
structure between males and females, we also calculated 
the codominant genotypic genetic distances by GenAIEx 
(Morin et al., 1994; Taylor et al., 1997) and the average 
relatedness (r) by COANCESTRY. If males are more 
related to proximate females, then inbreeding may occur 
(Lebigre et al., 2010). For the genetic distance, we expect 
the philopatric sex to have significantly higher genetic 
distance values than the dispersing sex. Additionally, we 
investigated the proportions of genetic distance with 
population B for different distance classes among males 
and females, to account for the sex-biased difference in 
spatial genetic structure. 


2.7. Spatial autocorrelation analysis Spatial genetic 
structure of P. shangchengensis was assessed using a 
multivariate spatial autocorrelation method, which 
strengthens the spatial signal by reducing stochastic 
noise resulting from the traditional allele-by-allele and 
locus-by-locus analyses (Smouse and Peakall, 1999). 
We calculated the squared genetic distance between 
two individuals following the method of Peakall and 
Beattie (Peakall and Beattie, 1995), and the sampling 
distance of each individual constituted the pairwise 
geographical distance. The autocorrelation coefficient (R) 
was calculated for a specified number of distance classes. 
The 95% confidence intervals for the estimate of R were 
generated by bootstrapping (1000 repeats). In this study, 
significant autocorrelation was declared only when 
r exceeded the permutation 95% confidence intervals 
around the null hypothesis of zero, and the bootstrap 
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95% confidence intervals around r did not exceed zero. 
Significant positive R implies that individuals within a 
particular distance class are genetically more similar 
than expected by random. Additionally, we conducted 
the Heterogeneity test to check the significant support 
of correlogram (Banks and Peakall, 2012). The spatial 
autocorrelation analysis was conducted including three 
parts (the entire sample with identified sexy males, and 
females) in GenAlEx. 


3. Result 


3.1. Basic population genetic information A total of 
195 adults of P. shangchengensis were obtained from three 
sampling sections of river, including 20 individuals for 
Population A, 153 individuals for Population B, and 22 
for Population C (Figure 2 and Table 2). According to sex 
characteristics, the number of identified individuals is 
187, with 101 male individuals and 86 female individuals 
(A: 10 male, 10 female; B: 79 male, 68 female; C: 12 male, 8 
female). 

In this study, small percentage of the microsatellite 
loci significantly deviated from HWE (Population A: 
PSH130, PSH463, PSH131; Population C: PSH133, 
PSH118) (Table 1). The frequencies of null alleles varied 


Table 2 Genetic diversity from three section of river (A, B and C 
population) in P. shangchengensis in China. 


n Na Hy + SE H; + SE Fs + SE 
A 32 10.17 OFT 023) 0.79 + 0.14 0.128 
B 156 15.50 0.71 40.19 0.83 + 0.16 0.145 
Շ 22 9.67 0.68 + 0.27 0.78 ੨ 0.15 0.151 


Notes: Sampling size (n), Number of alleles (Na), Observed Heterozygosity 
(Ho), Expected Heterozygosity (7), inbreeding coefficient (Fi). 


from 0.0094 to 0.26 and no loci showed evidence of 
null alleles. All loci showed wide variation in allele size 
range (Table 1). The numbers of alleles per locus (Na) 
ranged from 6 to 27 with the average of 15.83 (Table 1), 
indicating relatively high allelic diversity. The observed 
and expected heterozygosities (Ho, Hx) were 0.29-0.90 
(average 0.71) and 0.37-0.92 (average 0.83), respectively 
(Table 1). Moreover, the PIC values ranged from 0.33 
to 0.92 (average 0.81) (Table 1), indicating that these 
markers are highly informative. In total, 12 microsatellite 
loci showed polymorphism (Table 1). Among three 
populations, the number alleles across loci was higher in 
population B (15.50) followed by population C (11.33) and 
population A (10.17) (Table 2). Additionally, the higher 
values of observed and expected heterozygosities were 


| Vol. 11 


found in population B (Ho: 0.71; քլ: 0.83) and population 
A (Ho: 0.71; է: 0.79) followed by population C (Ho: 0.68; 
H,: 0.78). The replicated genetic diversity analysis for 
relatively equal sample sizes across the three populations 
(pop A: 20, pop B: 20, pop C: 22) based on 18 replicates 
showed very similar results, indicating that Population B 
had the higher genetic diversity, no matter Na or and Ho 
and H; (Table S1). 


The AMOVA 
analysis revealed a low but significant level of genetic 


3.2. Population genetic structure 


differentiation among the three populations (Fyr = 0.009, 
P = 0.001). Among the three populations, Population 
A was significantly statistical different from the other 
two populations (vs. population B Fy; = 0.010, P = 0.001; 
vs. population C 7੧ = 0.014, P = 0.002). Moreover, 
Populations B and C were also significantly statistical 
different (Fy; = 0.008, P = 0.003). The replicated AMOVA 
analysis at relatively equal sample sizes across the three 
populations (pop A: 20, pop B: 20, pop C: 22) based on 
18 replicates showed similar results, indicating little or 
no bias in estimation due to sample size variation (Table 
$2). The STRUCTURE results showed the significant 
genetic differentiation between population C and the 
other two populations (A and B) due to the proportion 
of genetic clusters and maximum AK (Figure 3 and 
Figure S1). On the other hand, when K=8, the genetic 
clustering composition in population A also has a certain 
degree of differentiation with the genetic clustering in 
population B, suggesting that there is also a significant 
genetic difference between these two populations (Figure 
3). FCA also showed that the individuals from three 
populations experienced significant separation trends 
(Figure 4), indicated substantial genetic differentiation 
among these three populations. In addition, this genetic 
differentiation pattern also got the support from the 
results of genetic distance (Pair-wise Jost's Dest and F’,.) 
between populations (Table 3). 

The relatedness estimates (r) and inbreeding coefficient 
(Fis) exhibit a linear relationship (Figure 5). The inbreeding 
coefficient values within populations were larger than that 
those between populations, suggesting close relatedness 
within populations (Figure 5). According to relatedness 
estimates (r) and the inbreeding coefficient 01) within 
the population, it seems that small populations will have 
larger inbreeding pressure. As for the three populations, 
Population B has a relatively small inbreeding coefficient, 
which indicates that individuals within the Population 
B might experience relatively longer dispersal distance 
(Table 4, Figure 5). In addition, the relatedness estimates 
(r) and the inbreeding coefficient (੧) between Population 
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Table 3 Pair-wise Jost's Dest (upper diagonal) and F'st (below diagonal) 
among three geographic populations of P. shangchengensis. 


A B Շ 
A 0.040 (0.004) 0.047 (0.032) 
B 0.004 (0.004) * 0.031 (0.018) 
6 0.005 (0.032) 0.003 (0.017) 


Note: Significant tests are indicated with asterisks (P Հ 005, ”P < 0.01); 
Inside the brackets is the P value. 


A and B were smaller than the equivalents between B and 
C, consistent with the than bigger genetic distance (F¢;) 
and population differentiation between population A and 
B (Table 4, Figure 5). 


Due to the limited 


number of identified sexed individuals in Populations A 


3.3. Dispersal pattern analysis 


(n = 14) and C (n = 12), all dispersal pattern analyses were 
conducted on Population B (n = 120). In our study, the 
Assignment Index correction (AIC) across all populations 
showed that female-biased dispersal (mean male AIC = 
0.026, SE = 0.198; mean female AIC = -0.031, SE = 0.225) 
(Figure 6). In addition, the mean of pairwise relatedness 
values between male individuals (mean r = -0.0405, SD = 
0.1115) was higher than that of female individuals (mean 
r --0.0404, SD = 0.1131), they both deviated significantly 
from zero, which further supported the female-biased 
dispersal pattern. The mean codominant genotypic 
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genetic distances between male individuals (22.701) was 
higher than that of female individuals (22.236), which 
also supported the female-biased dispersal pattern due 
to significantly higher genetic distance values of the 
philopatric sex compared with the dispersing sex. 

Spatial genetic autocorrelation analysis showed that 
the first distance class (0-70 m), the R value among males 
(R = 0.014, P = 0.111; Figure 7A) was higher than that 
among females (R = -0.0145, P = 0.929; Figure 7B) and 
among all individuals in B population (R = 0, P = 0.468; 
Figure 7C), which suggested that the autocorrelation 
value of the first distance class was not significant, but the 
autocorrelation value indeed was much larger than that 
in Figure 7B. Additionally, R value of the four distance 
class among males also show similar pattern. On the other 
hand, the Heterogeneity Test show that the correlogram 
of R value among males (Omega = 32.143, P = 0.006) 
was significant supported, but the other two show 
different results (females: Omega = 17.590, P = 0.129; all 
individuals: Omega = 19.163, P = 0.095). Combined these 
clues, the spatial autocorrelation analysis also hint the 
female-biased dispersal. 

Further examination of spatial distribution of 
codominant genotypic genetic distances with different 
genetic classes also showed that in the first distance class 
(0-50 m), the proportion of large genetic distance in 


female pairs (e.g. 3 or 4) was higher than male pairs and 


Շ 


Figure 3 Bayesian population structure analysis օք P. shangchengensis, on the basis of microsatellite genotyping data estimated by STRUCTURE. (A) 
output for K = 2, (B) output for K = 3, (C) output for K = 8. Different colors represent the proportion of each individual belonging to different putative 


clusters. 
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male-female pairs (Figure 8). Additionally, during those 


three distance classes (1.€., 301-350, 351-400, 401-451 3 ն լ z S 
m), the high genetic distances (e.g. 3 or 4 ) in male pairs 5 2 "Ը 
occupy a extremely low proportion; however, in the £ ի + 
female, it is the opposite situation. Theose results also 25 

hinted the female-biased dispersal in P. shanchengensis ի z 
(Figure 8). In addition, according to the relationship 111255 ; 18000 A 
percentage in population B, the percentage of female pairs 6000 -4000 -2000 200 4900 ge a A 
with genetic relationship was smaller than that of male min 10000 ਦਾਨਾ Է» 


Axis 1 (56.02 %) 
pairs, which also supported the female-biased dispersal 


(Table 5). In addition, following widely used criteria (0.25 


Figure 4 The factorial correspondence analysis across popula- 


Հ r < 0.5, moderate inbreeding; r > 0.5, high inbreeding) tions of P. shangchengensis. Three populations are indicated by 
(Hansson et al., 2007), we can infer that the most of different colors: red, population A; blue, population B; black, pop- 
ulation C. 


relatedness values within population B were smaller than 
0.25, which indicated that most of the individuals were 
out of inbreeding (Table S3). that in P. shangchengensis, smaller populations harbor 
lower genetic diversity and suffer higher inbreeding 
4, Discussion risk as one would predict. Studies have shown that 
negative consequences, such as low genetic diversity and 


4.1. Genetic diversity in isolated populations In inbreeding, are caused by increasing difficulty for species 


conservation biology, population isolation often leads to to access resources, such as food, mates, home domain 


a series of negative consequences, such as losing genetic size (Fahrig, 2003; Balkenhol et al., 2015). In this study, the 
diversity, reducing individual fitness, and increasing landscape barriers (the waterfall and underground river) 
the risk of inbreeding, which are closely related to the 
species’ viability and extinction risk (Ellstrand and Elam, 
1993; Avise and Hamrick, 1996; Fang et al., 2002; Keller 
and Waller, 2002; Burns et al., 2004; Reed, 2005; Collier 
et al., 2010; Hu et al., 2010). In this study, Populations 4.2. Landscape barriers and spatial genetic structure 


A, B and C are isolated from each other by a waterfall Due to limited dispersal or restricted gene flow, multiple 


in the stream habitat produced discontinuous sections 
and have led to the heterogeneity in the distribution of 


resources (e.g. food, mate, home domain). 


and underground river, Population A and Ը is smaller factors may influence spatial genetic structure, including 
than Population B. Overall, Population B showed higher geographic barriers, risk associated with dispersal, 
genetic diversity and a lower inbreeding coefficient than environmental heterogeneity, population density, 
Populations A and C (Tables 2, 4 and Table S1). It appears etc. (Koenig et al., 1992; Bowler and Benton, 2005; 


Table 4 The Statistical parameters of relatedness estimates (r) and the inbreeding coefficient (Fis) of each pair of individuals within and between 
populations in P. shangchengensis. 


A-A B-B C-C A-B A-C B-C 
r Fig r Fis r Fs r Fis r Fis r Fs 
Average 0.0017 0.0221 0.0341 0.0094 -0.0483 -0.0326 0.0426 0.0116 -0.0680 -0.0345 0.1772 0.0752 
SD 0.1346 0.0307 0.1757 0.1041 0.2454 0.1527 0.0525 0.0165 0.0659 0.0858 0.0523 0.0107 
Median -0.0029 -0.0067 0.0341 0.0094 -0.0483 -0.0326 0.0426 0.0116 -0.0680 -0.0345 0.1772 0.0752 
95% HPD Lower -0.0197 0.0172 0.0309 0.0075 -0.0799 -0.0523 0.0407 0.0110 -0.0743 -0.0427 0.1754 0.0748 
95% HPD Upper 0.0231 0.0270 0.0372 0.0113 -0.0166 -0.0129 0.0445 0.0121 -0.0617 -0.0262 0.1790 0.0756 


Table 5 The estimated relationship results by the highest likelihood within population B. Note: Ս, Unrelated; HS, Half Sibs; FS, Full Sibs; PO, Parent / 
Offspring. 


Relationship Proportion 
U FS HS PO Total U FS HS PO 
Male pairs 2807 9 265 0 3081 91.107 0) 292 8.601 0. 000 


Female pairs 2083 5 190 0 2278 91.440 0.219 8.341 0. 000 
Male-Female pairs 4925 11 435 1 5372 91.679 0. 205 8.098 0.019 
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Saastamoinen et al., 2018). Streams possess a dendritic 
hierarchical network structure, where smaller stream 
channels join to form larger ones in a dendritic pattern 
(Horton, 1945). For many stream-dwelling organisms, 
the consistent network architecture may form natural 


landscape barriers and may lead to local divergence of 


A Relatedness scatter 
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populations (Lowe et al., 2006; Campbell Grant et al., 2007; 
Campbell Grant et al., 2010) that ultimately face different 


evolutionary and ecological processes in different stream 


sections (Finn et al., 2006; Muneepeerakul et al., 2008; 
Grant et al., 2009; Mullen et al., 2010). For example, the 


Idaho giant salamander (Dicamptodon aterrimus) shows 


B Relatedness scatter 
y=0.5045x + 0.01 


R? = 0.4553 


D Relatedness scatter 
y = 0.4715 + 0.0052 
R? = 0.2784 


4 2= (0.216 
Fis 3 + 
2 
1 
-04 i 4 0.6 
r 
C Relatedness scatter 
y= 0.4507x + 0.0066 
Fis1.5 R? = 0.2493 
+ 1 
0.5 
0.6 
r 
E Relatedness scatter 
y = 0.4604x - 0.0027 
0.4 R = 0.5397 
Fis 


F Relatedness scatter 
y = 0.4974x + 0.0077 
R? = 0.4348 


Figure 5 The relatedness scatter of relatedness estimates (r) and the inbreeding coefficient (Իլ) within or among population in P. shangchengensis. (A) 
within population A, (B) within population B, (C) within population C, (D) between population A and B, (E) between population A and C, and (F) be- 
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Figure 6 Frequency distribution of corrected assignment indices (AIC) of P. shangchengensis. 


Asian 
Herpetological 
Research 


214 


70 140 210 280 350 420 
Distance Class (End Point, m; Male pairs) 


70 140 210 280 350 420 


70 140 210 280 350 420 


Distance Class (End Point, m; All ind pairs) 


Figure 7 Spatial genetic autocorrelation analyses among (A) male pairs 
(n = 3081 pairs), (B) female pairs (n = 2,278 pairs), and (C) all individual 
pairs (n = 10,731 pairs) in Population B. The genetic autocorrelation co- 
efficient (R) is a function of geographic distance, with the permuted 95% 
confidence intervals indicating random spatial genetic structure (dashed 
lines) and the bootstrapped 95% confidence error bars around r. U, upper 
95% confidence interval; L, lower 95% confidence interval. The smallest 
distance class size was set as 70 m. 


genetic differentiation among streams within catchments, 
consistent with a stream hierarchy model of population 
structure (Mullen et al., 2010). On the scale of this stream 
network, the patterns of genetic differentiation were 
often driven by the limited migration across the network 
architecture (Argentina et al., 2018; Jaisuk and Senanan, 
2018). 

At a smaller spatial scale, it is necessary to consider 
the influence of landscape factors on genetic pattern, 
especially landscape barriers that lead to discontinuous 
habitats, such as waterfalls and underground rivers in 
streams. Actually, the contiguous stream population that 
we assumed was composed of many smaller populations 
of different sizes (Jaisuk and Senanan, 2018). Whether 
landscape barriers substantially affect genetic structure 
needs further study. In this study, our data suggested 
that even at a smaller spatial scale, the main stream 
landscape barriers (waterfall and underground river) 
can still produce a clear signal of population genetic 
differentiation in P. shangchengensis (Table 3, Figures 3 and 
4). Populations A, B and C, are isolated by the waterfall 
and underground river. Population connectivity is critical 
to many ecological processes and conservation goals 
(Curtis and Bevers, 2002; Hayden, 2006; Saastamoinen 
et al., 2018). The key factor to population connectivity 
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is dispersal, which greatly affected population genetic 
structure, population dynamics, population persistence 
and adaptability (Clobert et al., 2001; Bowler and Benton, 
2005; Orsini et al., 2013; Smith et al., 2016; Olah et al., 
2017). Generally, limited by the characteristics of the 
linear habitat, the stream-type organisms are usually 
sensitive to the landscape barrier that effectively forms 
discontinuous stream habitat (Lowe et al., 2006; Campbell 
Grant et al., 2007; Campbell Grant et al., 2010). In this 
study, the genetic differentiation among Populations A, B, 
C indicate that waterfalls and underground rivers act as 
effective natural landscape barriers for P. shangchengensis. 
As a narrowly distributed species, our results suggest 
that the distribution of P. shangchengensis was affected by 
landscape factors, which may facilitate to understand the 
population divergence and evolutionary history of it. 


4.3. Female-biased dispersal Sex-biased dispersal may 
occur when the evolutionary pressures between sexes 
are different (Bonte et al., 2011). In recent years, sex- 
biased dispersal studies have been carried out in various 
vertebrate animal groups, such as fish (Bekkevold et al., 
2004; Stiver et al., 2007), amphibians (Palo et al., 2004; 
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Figure 8 Spatial distribution of codominant genotypic genetic distances 
among (A) male pairs (n = 3081 pairs), (B) female pairs (n = 2278 pairs), 
and (C) male-female pairs (n = 5372 pairs) in Population B. The distance 
classes was set as 50 m. The genetic distance between individuals (0, 1, 2, 
3 and 4) was calculated based on the method from Smouse and Peakall, 
(1999). 
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Liebgold et al., 2011; Wang et al., 2012), reptiles (Dubey et 
al., 2008), birds (Paris et al., 2016) and mammals (Héner 
et al., 2008). Among those groups, mammals usually show 
male biased dispersal and birds typically show female- 
biased dispersal (Greenwood, 1980; Pusey, 1987; Handley 
and Perrin, 2007). However, the dispersal patterns and 
mating strategies were diverse in amphibians and reptiles. 
For example, the direct-developing salamander Plethodon 
cinereus showed females as philopatric and male-biased 
dispersal due to the local mating competition (Liebgold et 
al., 2011); in the socially monogamous lizards of the genus 
Egernia, similar dispersal patterns were found (Gardner et 
al., 2001; Stow et al., 2001; Chapple and Keogh, 2005); In 
contrast, the dispersal pattern of the common frog (Rana 
temporaria) is female biased (Palo et al., 2004). Female- 
biased dispersal has also been observed in the Columbia 
spotted frog (R. luteiventris) (Pilliod et al., 2002) and the 
bullfrog (R. catesbeiana) (Austin et al., 2003). 

In this study, based on multiple analytical approaches 
(Figures 6, 7 and 8), including the Assignment Index 
correction (Alc), genetic distance, pairwise relatedness 
values of different gender (R) and spatial autocorrelation 
analysis, female-biased dispersal and male philopatry 
of P. shangchengensis was concluded. In nature, many 
animal groups show female-biased dispersal, such as 
fish (Taylor et al., 2003), amphibians (Pilliod et al., 2002; 
Austin et al., 2003; Palo et al., 2004), birds (Greenwood, 
1980; Paris et al., 2016; Pusey, 1987) and mammals 
(Hammond et al., 2006; Nagy et al., 2013; Tucker et 
al., 2017). In anurans and birds, the general trend for 
female-biased dispersal had been proposed (Greenwood, 
1980; Pusey, 1987; Pilliod et al., 2002; Austin et al., 2003; 
Palo et al., 2004; Wang et al., 2012; Paris et al., 2016). 
Unlike the amphibian species above mentioned, P. 
shangchengensis is a stream=salamander. In this study, 
we found differences in dispersal distances between 
the sexes, indicating female-biased dispersal and male 
philopatry of P. shangchengensis. Our results enrich the 
basic dispersal pattern of amphibian research, especially 
for the aquatic stream salamanders. Future studies 
should focus on whether female-biased dispersal acts 
as the general trend for aquatic Caudata. Unlike the 
female-biased dispersal of anurans (Austin et al., 2003; 
Lampert et al., 2003; Palo et al., 2004; Wang et al., 2012), 
it seems that there are no consistent sex-biased dispersal 
patterns in Caudata, although previous studies reported 
the male-biased dispersal pattern of two terrestrial 
salamanders (alpine salamander, Salamandra atra; red- 
backed salamander, P. cinereus) (Liebgold et al., 2011; 
Helfer et al., 2012). For Caudata, different habitat types 
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(e.g. terrestrial-type, stream-type) might have impacts on 
its dispersal patterns. Combined with the female-biased 
dispersal of P. shangchengensis, we think that there are 
multiple dispersal patterns (male-biased dispersal, female 
biased dispersal) corresponding to different habitat types. 
Therefore, for verifying the hypothesis, further research 
on dispersal patterns associated with habitat type may be 
needed. Additionally, the potential forces that facilitate 
or limit dispersal should also be studied to refine the 
understanding of the evolution of sex-biased dispersal. 
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Table S1 Genetic diversity across loci of population after shuffle in P. shangchengensis from three section of river (A, B and C population) in China. 
Pop Na Ho Hy Fs Pop Na Ho Hi Fs 


n 
էմ 
= 
Ծ 
Ի 
=] 
Յ 
Հ 
o 
© 
ԻՎ 
[e] 
> 
a 


n 
էյ 
= 
So 
a 
f=] 
[=] 
rox 
o 
© 
a 
f=} 
Ծ 
X 


an 
m 
= 
Ծ 
= 
5 
=] 
օօ 
=l 
ԷԶ 
ਕੰ 
Լյ 
Յ 
Հ 


a 
fest 
a 
o 
Ի 
5 
Յ 
Հ 
= 
= 
ਕੰ 
= 
> 
5 


n 
էյ 
=} 
Տ 
a 
So 
Տ 
ਹਾ 
o 
© 
- 
էթ) 
Ծ 
X 


n 
Էյ 
a 
© 
= 
Տ 
© 
- 
լա) 
է) 
A 
© 
੨ 
N 


2 
ԷՋ 
= 
o 
ps 
> 
Յ 
Հ 
=] 
= 
ਕੰ 
թ 
> 
a 


n 
ies! 
a 
= 
੦੦ 
So 
Տ 
ਹਾ 
=] 
է) 
K 
=] 
Ծ 
X 


n 
էղ 
an 
© 
=m 
լ) 
Տ 
- 
e] 
© 
Բ 
Տ 
Տ 
N 


n 
esl 
<< 
o 
& 
5 
Յ 
Հ 
o 
Ծ 
- 
f= 
> 
a 


੮3 
fest 


1.2 


լ) 
D 
ԵՀ 
So 
© 
a 
5 
Ծ 
E 


an 
fest 
= 
> 
= 
=] 
Ծ 
੦੦ 
=] 
ԷԶ 
Ի 
Յ 
Հ 


n 
էմ 
= 
Ծ 
Ի 
5 
Յ 
Հ 
[e] 
© 
ԻՎ 
> 
a 


an 
էմ 
= 
> 
੦੦ 
թ) 
> 
a 
=] 
© 
X 
o 
Ծ 
a 


੮3 
fest 
= 
Ծ 
= 
=] 
Ծ 
੦੦ 
੮ 
ԷԶ 
Ի 
=] 
Յ 
Հ 


n 
fest 
੩ 
Ծ 
bh 
=] 
Յ 
Հ 
=] 
ԷԶ 
Ի 
o 
> 
a 


n 
էյ 
_ 
= 
տ 
լ) 
է) 
X 
© 
© 
a 
Տ 
է) 
X 


n 
էղ 
<< 
© 
= 
լ) 
Տ 
Եշ) 
Տ 
© 
A 
© 
Տ 
N 


2 
fest 
= 
= 
A 
f=] 
Յ 
Հ 
=] 
© 
ਕੰ 
թ 
Տ 
a 


n 
fest 
Է 
a 
lon 
= 
=] 
a 
=] 
> 
mn 
=] 
Ծ 
տ 


੮3 
m 
੩ 
੨ 
= 
= 
=] 
օօ 
=] 
o 
A 
= 
Յ 
Հ 


an 
esl 
= 
Ծ 
Ի 
5 
੨ 
Հ 
o 
© 
ԻՎ 
o 
> 
a 


n 
ies 
= 
੨ 
Տ 
o 
Տ 
ਹਾ 
o 
© 
K 
o 
Ծ 
X 


n 
m 
(੪ 
੨ 
= 
5 
=] 
੦੦ 
=] 
ԷԶ 
Ք:Պ 
° 
Յ 
Հ 


n 
fest 
an 
Տ 
ՒՊ 
=] 
Յ 
Հ 
=] 
= 
ps 
= 
> 
a 


n 
էյ 
<< 
~ 
੧੩ 
Oo 
Տ 
ਹਾ 
o 
© 
a 
f=} 
Ծ 
X 


n 
੮ 
a 
© 
= 
Տ 
© 
- 
լա) 
է) 
A 
լ) 
© 
N 


n 
52) 
= 
© 
- 
Տ 
੨ 
Հ 
So 
o 
` 
So 
o 
an 


n 
m 
= 
Ջ 
Հ 
Տ 
[e] 
a 
Տ 
Տ 
տ 
լ) 
Տ 
տ 


n 
m 
a 
=] 
= 
Տ 
Տ 
co 
Տ 
է) 
Ի 
Տ 
fo} 
Հ 


n 
m 
= 
եռ) 
& 
Տ 
Տ 
Հ 
Տ 
է) 
Ի 
Տ 
© 
੦. 


n 
ies) 
Տ 
Տ 
੦੦ 
Տ 
੨ 
Հ 
So 
ծ 
` 
So 
Փ 
a 


n 
m 
= 
© 
= 
Տ 
Ջ 
੦੨ 
Տ 
ե) 
A 
լ) 
Ջ 
Հ 


n 
m 
= 
եռ) 
Ի 
Տ 
Տ 
N 
Տ 
է) 
Ի 
լ) 
Փ 
an 


n 
m 
Տ 
ԻԷ») 
w 
Տ 
Տ 
տ 
Տ 
է) 
Ի 
Տ 
ե») 
Ի 


n 
m 
— 
ծ 
= 
Տ 
© 
੦੨ 
Տ 
ե») 
Ի 
© 
Ջ 
Հ 


n 
եյ 
- 
© 
- 
Տ 
੨ 
Հ 
So 
o 
- 
So 
o 
a 


n 
m 
= 
m 
Ww 
Տ 
Տ 
Հ 
So 
Տ 
an 
[e] 
Տ 
Հ 


n 
m 
m= 
© 
= 
Տ 
Տ 
co 
Տ 
է) 
ਅੰ 
Տ 
Տ 
Հ 


n 
ies) 
- 
=] 
- 
Տ 
Տ 
Հ 
So 
o 
& 
So 
o 
a 


n 
m 
- 
Է. 
Տ 
© 
a 
լ) 
ե) 
X 
e] 
© 
տ 


n 
m 
pi 
ծ 
= 
Տ 
© 
co 
Տ 
է) 
& 
Տ 
fo} 
Հ 


n 
m 
= 
© 
Ի 
Տ 
jo} 
Հ 
Տ 
է) 
Ի 
Տ 
© 
੦. 


n 
m 
- 
N 
ԷՇ) 
Տ 
© 
a 
So 
© 
տ 
լ) 
լ) 
E 


n 
m 
— 
> 
= 
Տ 
© 
co 
Տ 
ե») 
Ի 
D 
Ջ 
Հ 


n 
m 
= 
եռ) 
Ի 
Տ 
Տ 
Հ 
Տ 
է) 
Ի 
Տ 
fo) 
© 


n 
m 
- 
Տ 
S 
a 
oS 
Տ 
տ 
oD 
© 
տ 


n 
m 
= 
> 
= 
Տ 
Ջ 
੦੨ 
Տ 
ե») 
Ի 
Տ 
Ջ 
Հ 


n 
ies) 
- 
=] 
- 
Տ 
Տ 
Հ 
So 
o 
- 
o 
Փ 
ԵԶ 


n 
m 
լա) 
ԷՋ) 
a 
Տ 
© 
տ 
լ) 
Տ 
տ 
oD 
© 
wo 


n 
m 
= 
oO 
= 
Տ 
Տ 
co 
Տ 
է) 
Ի 
Տ 
fo} 
Հ 


n 
ies) 
- 
© 
- 
լա) 
੨ 
Հ 
So 
ծ 
- 
So 
o 
ԵԶ 


n 
m 
- 
Տ 
© 
տ 
D 
Տ 
X 
լ) 
© 
տ 


n 
m 
= 
Ծ 
- 
Տ 
= 
co 
5 
o 
A 
=] 
੩ 
Հ 


Table S2 Results of Fy; after shuffle calculated from hierarchical AMOVA based on microsatellites of P. shangchengensis. 
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Figure S1 Plots of AK as a function of the number of putative 


groupings (K) in Bayesian population structure analysis of P. 
shangchengensis. 


